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Abstract. Using ideal and viscous hydrodynamics, the ratio of azimuthal moments V4,/{v2) 2 is 
calculated for pions, protons, and kaons in y / s=200 A*GeV Au+Au collisions. For any value of 
viscosity here is little dependence on particle species. Ideal hydrodynamics and data show a flat 
curve as a function of pt- Adding viscosity in the standard way destroys this flatness. However, 
it can be restored by replacing the standard quadratic ansatz for Sf (the viscous correction to 
the distribution function at freeze-out) with a weaker momentum dependence. 



1. Introduction 

The azimuthal distribution of detected particles in a relativistic heavy ion collision contains 
valuable information about the medium created in such a collision. This angular distribution is 
typically characterized by its Fourier components 
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where cj) is the angle of the outgoing particle momentum with respect to the collision plane for 
each event. 

The first non- vanishing coefficient is the elliptic flow coefficient V2, which is an important 
observable that is sensitive to bulk properties of the collision fireball and has been the object 
of much study jlj. Much less studied is the next coefficient, V4. In these proceedings, results 
are presented from ideal and viscous hydrodynamic simulations of Au+Au collisions at the 
Relativistic Heavy Ion Collider (RHIC). The main results of these calculations can be found 
in Ref. [2j, to which this can be considered a companion article. The details of the viscous 
hydrodynamic model can be found in Ref. [3]. 



2. Analytic prediction 

A simple argument by Borghini and Ollitrault [4J indicates that V4 is expected, at large pt, to be 
largely determined by V2- It is then most useful to construct the ratio ^4/(^2) 2 and investigate 
deviations from this "ideal hydrodynamic" scaling. 

To summarize, one can define a function u max ((f)), which is the magnitude of maximum 
fluid velocity on the freeze-out surface which is flowing in the azimuthal direction cj). For 



symmetric collision systems, this fluid velocity on the freeze-out surface will have only even 
Fourier components: 

u m a X (<^) = 17(1 + 2V 2 cos(2<£) + 2V 4 cos(4</>) + ...)• (2) 

If one plugs this expression into the ideal (zero viscosity) Cooper-Frye formula (neglecting 
quantum statistics) 

J 3 N f ,„„ f p^u^ 



performs a saddle-point integration, and expands to leading order in these "intrinsic" elliptic 
and quadrangular flow coefficients V 2 <C 1 and V4 ^ 1, the result is 

nipt) ( pt ~ mtV "> ^ 
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with m t = ypt+m? and v = 17/ \f\ + U 2 . This suggests that it is useful to construct the ratio 

Vi/{v 2 ) 2 : 

va __ 1 VaT = 1 V± 1 

(v 2 y 2 + V 2 U (p t - m t v) 2 + V 2 v 2 {p t ) ' 1 j 

So, according to this prediction, the ratio should approach 1/2 at large pt, with 1/pt corrections 
at smaller transverse momenta. Note that this prediction does not depend on any details of 
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Figure 1. Charged hadron Vi/(v 2 ) 2 as a function of transverse momentum, using Glauber 
and CGC initial conditions with impact parameter b = 10 fm, and using the standard quadratic 
ansatz for Sf. Data are from STAR (40-50% centrality) [S] and PHENIX (45-50% centrality) [§], 
scaled down by a factor of 1.38 to account for fluctuations 



the collision, as long as the assumptions are satisfied (that V4 and V2 are small and that one 
is considering large enough pt compared to the freeze out temperature). Also note that, when 
plotted as a function of V2 instead of pt, there should be no dependence on the particle species. 

3. Results 

Results from ideal hydrodynamics confirm these expectations. As detailed in Ref. [2], when the 
viscosity is set to zero, the ratio v^/^2) 2 approaches roughly 1/2 at large pt with a correction 
that behaves as 1/pt, as predicted from Eq. (|6|). Both the magnitude and sign of this correction 
depend strongly on the freeze-out temperature. A freeze-out temperature close to Tf = lAOMeV 
(the temperature that gives the best fit to other data [3]) gives a flat dependence on p t , which 
is what is seen experimentally (see, e.g., Fig. [1]). 

Using different choices of initial conditions for hydrodynamic evolution [Glauber-type or 
Color-Glass-Condensate (CGC)-type] has a small effect, indicating little sensitivity to the initial 
eccentricity. This is in contrast to v% itself, which is largely driven by the initial eccentricity. 

As predicted, and as seen experimentally 019], there is no dependence on particle species at 
high p t . Since there is little dependence on transverse momentum, this is true whether v^jiv^) 2 
is plotted versus pt or V2 ■ 

Data, however, show a value that is closer to 1, for a large range of centrality and pt- Most 
of the discrepancy can be understood by noting that t>4 and vi are obtained separately when 
averaged over events. Event-by-event fluctuations in V4 and V2, which are not included in these 
hydrodynamic calculations, then increase the ratio v^/(v2) 2 [I]. 

To add viscosity to the calculation, one must first make a choice for 5f, the viscous correction 
to the equilibrium distribution function at freeze-out. Until recently, all groups used a quadratic 
ansatz for the momentum dependence of Sf, but the correct form is not known, and it can have 
a significant effect |lUj . 

Adding viscosity and using the standard quadratic ansatz results in a significant change to 




Figure 2. Identified particle v^/^2) 2 versus V2 from viscous hydrodynamics with CGC initial 
conditions and rj/s = 0.16 at b = 8 fm using the standard quadratic ansatz as well as a linear 
ansatz for the viscous correction, 5f, to the equilibrium distribution function at freeze-out [10]. 



the pt dependence of v^/{v2) 2 , as shown in Fig. [TJ There is again little dependence on the 
initial conditions. Although there is some dependence of the overall size of v^jiv-^) 1 on initial 
conditions, the shape of the curve is largely determined by only the viscosity. In addition, 
a large centrality dependence is introduced (not shown). Figure [2] shows that, as with ideal 
hydrodynamics results, there is little dependence on particle species at high p t when plotted as 
a function of V2, even when the viscous corrections are large, and regardless of the choice for Sf. 

Using a weaker momentum dependence for 8f can make the pt dependence of ^4/(^2) 2 more 
flat (see Fig. [2D, and at the same time reduce the dependence on impact parameter. This may 
indicate that the standard quadratic ansatz is not correct and could give insight into the behavior 
of the hadron gas present at freeze-out (see Ref. [10J for more details). 

4. Conclusions 

Viscous hydrodynamic simulations of heavy ion collisions were performed, and the ratio of 
azimuthal moments v^/fa) 2 was computed. It was found that, as in ideal hydrodynamics, 
identified particles in viscous hydrodynamic simulations have the same value, except at low pt, 
when plotted as a function of v%. In addition, the shape as a function of pt depends little on the 
initial conditions, but depends strongly on the shear viscosity to entropy density ratio rj/s, when 
implemented using the standard quadratic ansatz for the viscous correction to the equilibrium 
distribution function at freeze-out. In contrast, a weaker momentum dependence for 5f can 
restore the flat dependence on transverse momentum that is seen in ideal hydrodynamics and 
data, and reduce the strong centrality dependence, which is not seen in experiment. 
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